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■ We study the problem of folding of the regular triangular lattice in the presence of 

3 \ bending rigidity K and magnetic field h (conjugate to the local normal vectors to the 
triangles). A numerical study of the transfer matrix of the problem shows the existence 
^ ' of three first order transition lines in the (iiT, h) plane separating three phases: a folded 
phase, a phase frozen in the completely fiat configuration (with all normal vectors pointing 
up) and its mirror image (all normal vectors pointing down). At zero magnetic field, a first 
order folding transition is found at a positive value K^, ~ 0.11(1) of the bending rigidity, 
corresponding to a triple point in the phase diagram. 
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1. Introduction 

It is tempting to try to describe geometrical objects like polymers (ID) or membranes 
(2D) in analogy with spin systems. Natural spin variables are provided for instance by 
the local normal or tangent vectors to the object, while elastic properties like bending 
rigidity naturally translate into some nearest neighbor spin coupling. However, the cor- 
respondence between geometrical objects and spin systems can be subtle, especially in 
two dimensions, where geometric constraints on, say, the normals to the membrane imply 
local constraints on the associated spin variables. Such constraints play a crucial role for 
tethered membranes, i.e. 2D polymerized networks with fixed connectivity, since they 
induce a crumpling transition by stabilizing an ordered phase in a region where the un- 
constrained spin system would be disordered. This phenomenon was recognized in [^-j^, 
where a continuous crumpling transition is predicted. 

Such a drastic change of statistical behavior is observed in the present paper, where we 
consider a spin system describing the thermodynamics of folding of the regular triangular 
lattice, a problem first considered in 0. Considered as a geometrical object, the lattice 
describes a tethered membrane skeleton, made of rigid bonds along which folds can be 
performed. Here we consider only complete foldings which result in two-dimensional folded 
configurations of the membrane. In such a process, each bond serves as a hinge between 
its two neighboring triangles, and is in either one of the two states: folded (with the 
two neighboring triangles face to face) or not (side by side)0. A folding configuration 
(folded state) of the system is entirely specified by the list of its folded bonds. This 
definition corresponds to a "phantom" membrane, where the folding process may imply 
self-intersections, and where one cannot distinghish in the folded state between different 
piling orders for superimposed triangles. 




-I- 



(a) (b) (c) (d) (e) (f) 




(g) 



(h) 



(i) 



(j) 



(k) 



Fig. 1: The eleven local fold environments for a vertex. Folds are represented 
by thick lines. One of the two possible spin configurations on the triangles is 
also indicated. 



^ We refer the reader to [^[Q for a more formal definition of folding. 
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With this simphfied definition, our folding problem can be formulated as an eleven 
vertex model, expressing that the immediate surroundings of a vertex in a folded state must 
be in one of the eleven local configurations depicted on Fig.l. In spite of its local definition, 
folding is a highly non-local operation. As explained in , we note that whenever a bond 
is folded, say, on the left half of a vertex, then another bond is folded on the right half, 
hence folds propagate throughout the lattice. 

In a previous work 0, we have computed the exact thermodynamic entropy per 
triangle which counts the number Zf of folded configurations for a finite lattice made of 
Nt triangles for large Nt —>■ oo. This was done by mapping the eleven vertex model above 
onto the three-coloring problem of the triangular lattice bonds, solved exactly by Baxter 
in its dual version, the three-coloring of the hexagonal lattice. The result reads 

s = lim ^Log Zf = Log q 

oo Q _ 1 /q (1-1) 

^3p(3p - 2) 2n ^ ^ ' 

As mentioned above, we will rather use here the alternative description of the model 
in terms of Ising spin variables a = ±1 defined on the triangles, and indicating whether 
they face up or down in the folded state. One can think of the spin as the normal vector to 
the triangle. Spin configurations are given together with the fold configurations on Fig.l. 
Note that there are two spin configurations for each folded state, due to the degeneracy 
under reversal of all spins, hence the partition function Z of the spin system is twice that 
of the eleven vertex model: Z = 2 Zf. It is clear from Fig.l that the only allowed vertex 
environments are those with exactly 0, 3 or 6 surrounding up spins. In order for a spin 
configuration to correspond to a folded state, the six spins cr^ around any vertex v must 
satisfy the local constraint 

Ey= ai = mod 3 , (1.2) 

i around v 

since = 2( number of up spins) — 6 is a multiple of 3 iff the number of up spins itself 
is a multiple of 3. Beyond the above counting of the number of allowed constrained 
spin configurations, it would be desirable to understand the effect of a bending energy 
for the folds, characterizing the rigidity of the membrane. In the spin language, this 
means the presence of a ferromagnetic Ising-like interaction energy —Jaiaj between nearest 
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neighbors. Most properties of the folded tethered membrane can in fact be investigated 
by studying the magnetic behavior of our constrained Ising spin system. The average 
magnetization M of the system is indeed an order parameter which is characteristic of the 
flatness of the membrane (|M| > for a configuration which is flat in average, and M = 
for a conflguration which is folded in average). This suggests to introduce a magnetic fleld 
H in the system (with energy —Hai per triangle), with no direct physical meaning for the 
membrane, but instrumental in revealing information on its average state of folding. In 
the following, we will therefore consider the constrained Ising model with Hamiltonian 

Ti-ising = -J'^(^i(^j - H'^ai . (1.3) 

For convenience we will use the reduced coupling and magnetic fleld 

K = J/keT ; h = H/keT . (1.4) 

In our study of this model, we will give numerical and theoretical evidence for the existence 
of a first order transition line in the {K, h) plane, between an ordered phase M = 1 (for 
h > 0) where the membrane is completely flat, and a disordered phase M = 0, where the 
membrane is folded and has a non-vanishing entropy. The most surprising fact is that we 
flnd no intermediary magnetization of the system in the thermodynamic limit. For h = 0, 
a flrst order folding transition still takes place, at a critical value Kc of the reduced Ising 
coupling K. All these results clearly show a drastic modiflcation of the thermodynamics 
of the standard Ising model, emphasizing the special role played by the constraint (|1.2| ). 

The paper is organized as follows. In section 2, we describe the transfer matrix that we 
shall use for numerical simulations on the thermodynamics of the constrained spin system, 
and show how to take advantage of some particular properties of this matrix. The results 
for the magnetization in the presence of a magnetic fleld are discussed in section 3, and 
lead us to formulate the abovementioned two phase (M = 0, 1) hypothesis. Under this 
assumption, we also derive a simple argument to calculate the critical value of the magnetic 
fleld where the transition between these phases takes place. Section 4 is dedicated to the 
precise study of the flrst order transition line in the thermodynamic limit. In particular, 
we flnd the critical value Kc beyond which the M = 1 phase persists even at zero magnetic 
fleld. We discuss the general phase diagram of the system in section 5, and gather more 
evidence for the flrst order character of the transition. Related topics are discussed in 
section 6, including the exact solution for the square lattice as well as some predictions of 
a possible antiferromagnetic transition within the M = phase, for negative K. Section 7 
is a brief conclusion. 
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2. Transfer matrix description 

We consider the folding of an infinite strip of triangular lattice of finite width L, with 
free boundary conditions on the edges of the strip. Imposing periodic boundary conditions 
at infinity, the partition function of the model is expressed as 

Z^^\K,h)= lim rTr(T(^)(K,/i)^)l^, (2.1) 

where T^^^ denotes the transfer matrix, acting as an operator on a column of size L, whose 
state is specified by the 2L spin values ai = ±1, z = 1, .., 2L on the triangles. 
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Fig. 2: The transfer matrix T*^"^) . 

The matrix element of T^^'^ between two consecutive columns, as depicted on Fig. 2, 
reads 

where T^^\0,0) imposes the local folding constraint (|1.2| ) on the six spins surrounding 
each of the L — 1 inner vertices (marked by black dots on Fig. 2 ), and and V^^^ are 
the usual temperature and magnetic field contributions to the transfer matrix of the Ising 
model, namely 

L-l 

^K},{<t}(0' 0) = n + + cr2^+2 + (J2,-i + a'^i + ^ai+i mod 3) 



1=1 



K 



2L-1 



i=l 
2L 



i=l 
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with S{x) the usual Kronecker delta function on integers. 

In the large limit, the partition function ( |2.1| ) is dominated by the largest eigenvalue 
of T^^\K,h), which we denote by X^IxiK, h), and the corresponding free energy per 
triangle reads 

= ^LogZW(i^,/i) = ^LogX^J:l^{K,h) . 

The thermodynamic free energy per triangle is then defined as the L —>■ oo limit 

-F{K,h)= \im ^Log Xi^liK, h) . (2.2) 

Our main task is therefore the numerical extraction of this eigenvalue, together with the 
corresponding eigenvector Vmax{K,h). We will also be interested in evaluating the sub- 
leading eigenvalue X^^l{K, h). 

The analysis can be simplified in view of the following properties of the transfer matrix. 
Due to the local constraint on the spin variables, the matrix T'^^\K, h), of size 2"^^ x 2^^, 
is sparse. The exact number Nl of non-zero entries of T'^^^ is easily computed by using 
the vertical transfer matrix of size 2, T(2)(0,0), namely 



t(2)(o,o)^ 



J {a'}, {a} 



Diagonalizing T*^^) exactly, we find 



Nl = - 
2 



1 /7 + v^\^ r 11 + 3v^l /7- v^\^ ril-3v^ 



17 J V 2 / I V17 



The first values of Nl are given in Table I below. A^^ grows like (5.56..)^ < 16^, which 
justifies the use of a diagonalization algorithm adapted to sparse matrices: starting from a 
given vector, we compute recursively its iterated image by T, which we moreover normalize 
at each step; the process converges to the eigenvector Vmax associated to Amax- The 
subleading eigenvalue Agub is obtained by use of the same algorithm, together with a 
suitable projection procedure, guaranteeing at each step that we substract the component 
along the maximum eigenvector. This projection procedure is made easy by the following 
symmetry property of the transfer matrix. Consider the picture of Fig. 2, and rotate it by 
180 degrees. One gets an identity between two transfer matrix elements, namely 



o'2L,---,o-2,a-i\(j^^,...,cj^,a-^ 
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This can be recast into T^^\K,hY = RT^^\K,h) R, where R is an involution, with 
matrix elements 



2L 



O-J^ ,(T2 , . . . jCTj^ ;(Tl ,CT2 , . . . ,0-21, 



1=1 



Any eigenvector v with eigenvalue A < Amax is then clearly orthogonal to Rvma.x- To 
project out the Vmax component off a given vector w, one simply has to perform the 
substitution 

(^^'maxk) 



W 



W 



{R 

From the knowledge of the eigenvector Vmax{K, h), we can get the expectation value of 
local observables of the system. For instance, the magnetization M is obtained as follows. 
The magnetization operator ^ acts diagonally on the columns of 2L spins, with diagonal 
elements A*{ct},{ct} = Si=i ^i- the N —>■ oo limit, its expectation value M = (//) reads 



{Rvi^l{K, h)\^iv'^L{K, h)) 



M = 



{Rv^rtUK,h)\v'r^L{K,h)) 



(2.3) 



As a test of the precision of our algorithms, we first compute Ainax(0,0) for various 
sizes L = 2, 3, .., 9. The results are summarized in Table I. 



L 


Nl 


A^ix(0,0) 


2 


88 


5.56155281 


3 


488 


7.89397081 


4 


2712 


11.32598261 


5 


15080 


16.34742307 


6 


83864 


23.67855022 


7 


466408 


34.37351897 


8 


2593944 


49.97256996 


9 


14426344 


72.72483625 



Table I: The number Nl of non-vanishing transfer matrix elements and 
the maximum eigenvalue Amax for strips of width L = 2, 3, .., 9. 
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This leads to an estimatell for the thermodynamic entropy per site — -F(0, 0) = Logg, with 
q = 1.208.., in very good agreement with a previous numerical estimate 0, and with the 
exact result (|1.1| ). 



3. Magnetization, critical magnetic field and the two phase hypothesis 

As mentioned above, the magnetization M of the system is an order parameter for 
the flatness of the lattice. At K = 0,h = the system is folded in average, with a 



non-vanishing folding entropy —F{0,0), therefore the magnetization M (|2.3|) vanishes 
identically. On the other hand, we will have M — > 1 for sufficiently large h > 0. 
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Fig. 3: Magnetization M versus magnetic field h for L = 4 and 6. The four 
curves correspond respectively (from the right to the left) to exp(i^r/2) = 
1 {K = 0) , 1.0333.., 1.0666.. and 1.1. The dashed vertical lines indicate the 

critical magnetic field hi^^ as predicted by eqn.(3.1). 



Fig. 3 represents the magnetization M versus the magnetic field h > 0, computed 
through the formula (|2.3| ) for strips of width L = 4 and 6, for several values of K, in the 
range [0, 0.2]. It clearly appears, and even more so for the larger L, that the magnetization 
tends to remain zero over a finite range of (small enough) magnetic fields, and then abruptly 
jumps to 1 when h reaches some critical value hi^"^ {K). Moreover, this critical field hi^"^ (K) 



^ We obtain the value of q as the limit of the sequence ql 
Aitken delta-2 algorithm (exponential fit). 



/Amax, extracted by the 
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is maximal for K = and decreases with increasing K. For a given K, hT\K) also 
decreases with increasing L, eventually reaching its thermodynamic limit hc{K) for L — 
oo. 

This somewhat unexpectedly rapid change in M is the first tangible sign of the ex- 
istence of a first order magnetic transition in the system. At this point, it is reasonable 
to infer that in the thermodynamic limit L —>■ oo, and for a given coupling K > 0, the 
magnetization M is exactly zero for < h < hc{K), and exactly one for h > hc{K), being 
therefore discontinuous at h = hc{K), with hc{K) a decreasing function of K. Indeed all 
our results in the following will corroborate this picture of a first order transition between 
two phases with respectively M — and M = 1, without any possible intermediate value 
of the magnetization. 

A first check of the above two phase hypothesis is actually provided by the possible 
derivation of the value of the critical magnetic field through the following simple theoretical 
argument. Suppose that for large enough but finite strip width L, we can already describe 
the system in terms of two phases M = and M = 1. In the phase M = {h < hi^\K)), 
the system is insensitive to the value of the magnetic field h, its partition function is 
therefore given by Z^^\K,h) ~ Z^^\K,0) = exp ( - 2LF^^\K,0)). In the fiat phase 
M = 1 {h > hi^\K)), the partition function is that of the pure state with all spins up, 
hence reads Z^^\K^ h) ^ exp ((3L — 1) iiT + 2L /i) , since 3L — 1 bonds separate the 2L 
triangles. The phase transition is then predicted to occur at the critical value of the field 
h where the free energies of the two phases are identical, namely 

/,W(K) = _^W(K,0)-|(3-i) , (3.1) 

where —F^^\K,Q) = ^Log X^a^iK , 0) can now be calculated numerically, directly from 
the zero magnetic field transfer matrix. The corresponding predicted values of the critical 
field are represented on Fig. 3 by dashed vertical lines, for the various values of K and L. 
The agreement with the observed transition point is excellent. In the following, we shall 
therefore consider h\^\K) of eqn.( |3.1| ) as giving the exact location of the transition point. 
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Fig. 4: The critical magnetic field hi^\K) resulting from eqn.( p.l|) , and the 



purported thermodynamic limit hc{K) (dashed line). 
4. Transition line and critical coupling K^, 

We are now interested in understanding the thermodynamic critical line hc{K) sepa- 
rating, in the {K, h) phase diagram, the M = and M = 1 phases. 

On Fig. 4, we have represented the curves hi^\K) as given by eqn. (|3.1|) for L ~ 
2,3, ...,8 and K G [0,0.36]. These curves enjoy the following properties. Eqn. (|3.lD ex- 
presses the critical field hi^\K) as the difference between the total zero field free energy 
and that of a particular state (the flat M = 1 state): hi^\K) is therefore positive. It 
is also clear that hi^\K) when i^T ^ oo as the system gets fully ordered for strong 
coupling. Finally hi^\K) is a decreasing function of K, indeed 

\K) 1 



dK 



1) 



(4.1) 



1 / In 



where the sum extends over the 3L — 1 active bonds (ij) of a column of 2L triangles, and 
E^^\K, 0) is the zero field average energy per triangle of a strip of width L. As {(JiCTj) is 
always smaller or equal to one, we deduce that hi^^ (K) decreases with increasing K. 
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When L — > oo, the curves hi^\K) of Fig. 4 tend to the thermodynamic critical hne 
hc{K). The properties mentioned above for finite L naturaUy extend to the thermodynamic 
Umit. The critical field, now given by 

K{K) = -F{K,Q)-^K , (4.2) 

is thus positive or zero. At -ftT = 0, /ic^'*(0) is nothing but the entropy per triangle 
-F(^)(0,0) = 2xLog A^4x(0,0), with A^ax(0,0) given in Table I. Consequently, in the 
thermodynamic limit L — > oo, we have the exact result 

/ic(0) =Logg= .189... , (4.3) 

the transition thus taking place at a non-zero value of h. By continuity, hc{K) will remain 
strictly positive for small positive K. Differentiating eqn. (|4.2| ), we get 

dhc{K) 3 

^^ = 2((^^-.)(^,o)-l) ^^^^^ 

where {aiaj) denotes the correlation function of any two neighbouring spins, and E{K, 0) 
is the zero field average energy per triangle of the system. Again, as {cfiCTj) is always 
smaller or equal to one, we deduce that hc{K) is a non-increasing function. Note that 
eqn. (|4.4|) is nothing but the Clapeyron relation 



dK ' ' Ml -Mo' 

where Ei (resp. Mi) denotes the Ising bending energy (resp. magnetization) in the phase 
z = or 1 on each side of the critical line. 

Once hc{K) is known, the h dependence of the system is determined since 

In turn from eqn.( ^.2| ), hc{K) is encoded in the zero field free energy —F{K, 0). 

A new interesting phenomenon can be read off Fig. 4. For large enough L, the decrease 
to zero of hi^\K) with increasing K takes place over a finite interval [0, i^c^^], with hi^^ ^ 
for K > Ki^\ This is best seen on Fig. 5 which represents the curves 
E^^\K,0)-{3L- l)/2LforL = 2, 3, 8 and G [0,0.36]. One clearly sees a jump in the 
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Fig. 5: dh[^\K)/dK = E^^\K,0) - (3L - 1)/2L, for L = 2,3, ...,8 and 
K G [0,0.36]. 

sfope of hi^\K) from a finite value (~ —1.2 for L = 8) to 0. Moreover the intersections 
of the various curves provide us with estimates of the critical values Kc^\ The latter 
decrease with L and converge to a limiting value Kc ^0.1. 

In the thermodynamic picture, this means the existence of a critical value of the 
coupling K, such that hc{Kc) = (i.e. — F(Kc, 0) = f-ftTc), and thus hc{K) = for all 
K > Kc- On Fig. 4, we have represented in dashed line the curve hc{K) as obtained by 
extrapolating! to L = oo the values of hi^\K) at fixed K. This direct extrapolation 
confirms the emergence of a critical Kc and predicts a value Kc = 0.11(1). From Fig. 5, it 
corresponds to E{Kc, 0) ~ 0, i.e. a vanishing nearest neighbor spin correlation (ciaj) ~ 0. 

^ Like for our numerical estimate of q, we obtain the value hc{K) as the logarithm of the limit 
of the sequence qL{K) = exp exp Lhc^\K), extracted by the Aitken delta-2 

algorithm. 
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Fig. 6: Phase diagram in the {K, h) plane. Three first order hnes h = 
hc{K), —hc{K) {K < Kc) and h = {K > Kc) separate the three phases 
M = 0, ±1 and meet at the triple point {K^ 0). The dashed line represents 
a constant magnetic field line, which crosses the transition line h = hc{K) at 
a critical value Kc{h). 

5. Phase diagram 

Our results are summarized on Fig. 6, representing the phase diagram of the sytem in 
the {K, h) plane. We extended the range of h and K to include real negative values. The 
phase diagram is clearly symmetric under h —>■ —h, while a negative K simply corresponds 
to an antiferromagnetic Ising coupling, which favorizes folding. Three first order transition 
lines h = hc{K), —hc{K) {K < Kc), and h = {K > Kc) separate the three phases 
M = 0, ±1. The line hc{K) naturally extends to negative K through eqn.( |3.1| ). For large 
negative K, the M = phase is dominated by the two pure antiferromagnetic states with 
alternating spins, with free energy —F{K, 0) hence hc{K) ~ —3K for K — > — oo. 

We expect nothing special to occur at i^T = for the line hc{K), where we have the exact 
result hc{0) = Logq = .189.... Instead our study predicts the existence of a triple point 
{Kc,0) at the positive value Kc — .11(1). For the physical zero magnetic field membrane 
problem, this corresponds to a first order folding transition at K = Kc- Within the domain 
M = 0, the system is insensitive to the magnetic field h. Along the constant magnetic 
field dashed line of Fig. 6, with K increasing from — oo, the free energy of the system is 
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—F{K,0) until one reaches the critical line K = Kf.{h) (inverse oi h = hc{K)), beyond 
which the free energy becomes the linear function —F{K, h) = + h. 




Fig. 7: Plot in dashed hnes of ji^ogX[^^^ sub(-^' ^) L = 6, exp{h) = 1.1 
(short dashes) and exp(/i) = 1.2 (long dashes), and K G [0,0.15]. The solid 
lines represent the free energies —F^^\k, h) = h + nK/12, and —F^\k, 0). 
The crossings occur at the critical couplings i^c^'*(Log 1.2) ~ 0.05 and 
K^^^ (Log 1.1) ~ 0.10. 

As usual for first order transitions in the transfer matrix formalism, this change of 
behavior results from the crossing of the two largest eigenvalues Amax and Asub of the 
transfer matrix T in the thermodynamic limit. Indeed this is already visible for finite L, as 
exemplified on Fig. 7, where we plot for L = 6 (in dashed lines) the two leading eigenvalues 
of T'^^\K,h) for two different positive values of h, in the vicinity of the corresponding 
critical points K^\h). One clearly sees the exchange of the two eigenvalues with 

K>K'^''\h) 



2L ^ 1 Asub 
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where -Ff>(K,h) = h + {3L - l)K/2L, and -Ff\K,h) = -F(^\K,0) denote respec- 
tively the free energy in the phase with magnetization M = 1 and 0, represented for L = 6 
on Fig.7 in sohd hnes. 



1.5 




-0.5 ' ' ' ' ' ' ' ^ 
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K 



Fig. 8: Finite (dashed) and zero (sohd) field energy versus Ising coupling, 
for L = 6 and exp(/i) = 1.2, 1.1, 1 respectively, and K G [0, 0.36]. 

The change of behavior of the free energy along a constant-field line at Kc{h) is 
confirmed at finite L by the plot of its derivative w.r.t. K, the Ising energy E'^^\K, h) = 
—OkF^^^K, h), which we represent on Fig.8 for exp(/i) = 1, 1.1, 1.2 and L = 6. We clearly 
see that the finite field energies E^^\K, h) exactly match the zero field energy E^^\K, 0) 
before K reaches the critical value K^\h), where they abruptly jump (the more so for 
higher h) to the asymptotic value (3x6 — l)/(2 x 6) = 17/12. 

The corresponding specific heat ci^\K,h) = dxE^^^K, h) is represented on Fig. 9 
for L = 6 in dashed (exp(/i) = 1.1, 1.2) and thick solid (exp(/i) = 1) lines. As expected, the 
finite field specific heats match the zero field specific heat until K reaches Kc^\h), where 
they exhibit a delta-function peak, before immediately (the more so for higher h) reaching 
a zero value. The peak for /j, = seems to be qualitatively different from that for finite 
h. To see how this peak develops as the size L grows, we have represented on the same 
figure the zero field specific heat for L = 2, 3, 8 in thin solid lines. The smoothness of 
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Fig. 9: Specific heat ci^\K, h) for exp(/i) = 1 (thick sohd fine), 1.1 and 1.2 
(dashed hnes) and ci^\K, 0) for L = 2, 3, 8, for K e [0, 0.36]. 

the curves might be the sign of some extra divergence on top of the delta-function in the 
thermodynamic specific heat Cv{K, 0), possibly of the form C^{K, 0) ~ {Kc — K)~" when 
K approaches the triple point from below. Our data do not allow us at present to reach 
any conclusion on this point. 

6. Discussion 

The main result of this paper is the phase diagram of Fig. 6. In particular, the ex- 
istence of a triple point {K = Kc > 0, /i = 0) at the boundary of the three phases with 
magnetization M = 0, ±1 is a non-trivial outcome of our study. It is interesting to note 
that a very similar picture can be obtained exactl^ in the case of the folding of the square 
lattice. As already noted in the thermodynamic entropy of folding — -Fsquare(0, 0) of 

An exact solution can be obtained for instance by using a transfer matrix which is diagonal 
in terms of folded line variables. 



i 



exp(li)= 
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the square lattice vanishes, due to the very strong constraint that, for this particular lat- 
tice, folds must propagate along straight lines, all the way through the lattice. A folded 
state of the square lattice is entirely specified by the data of its folded horizontal and 
vertical lines. For a square lattice of size L x L, this leads to a free energy per square 
— -^square (0, 0) = HniL^oo ;p-Log 4^ = 0, as annouced. Like in the triangular case, the 
square lattice folding problem is easily transformed into a (face) square lattice Ising spin 
system, with the local constraint that there are exactly 0, 2 or 4 spins up around each 
vertex. The zero field thermodynamic free energy —Fsq^^aTe{K,0) per square at an arbi- 
trary value of the reduced Ising coupling K is easily obtained as follows: at K = 0, it 
vanishes; for K — > oo, it tends to the fiat state value 2K; for K — oo, it tends to the 
completely folded state value —2K. From the usual convexity property of —F, we conclude 
that necessarily 

-Fsquare(i^,0) =2|K| . 

The free energy appears here simply as a competition between the contribution of the 
completely folded state -Fo(K, 0) = -2K and that of the flat state -Fi{K,0) = 2K, 
with — -Fsquare(-f^5 0) = Max(— Fq, —Fi). Aualogously, in the presence of a finite positive 
magnetic field h, the free energy results from the competition between the contribution of 
the completely folded state, insensitive to /i, —Fo{K, h) — —2K, and that of the flat (all 
spins up) state —Fi{K, h) = 2K + h. This yields 

-Fsquare(i^, h) = Max( - Fo(i^, h), -F,{K, h)) = I ^-^^^ ^ < . 

Consequently the system undergoes a first order phase transition from the completely 
folded state with magnetization M = to the fiat state with M = 1, along the critical line 

/ic,square(i^)=2(|i^|-i^) , 

the square lattice analogue of eqn. (|3.1|) . 

As shown on the resulting phase diagram of Fig. 10, the qualitative behavior of the 
system is very similar to the triangular case discussed above, corroborating a posteriori 
the diagram of Fig. 6. A crucial difference is that the triple point now sits at the origin of 
the {K, h) plane, i.e. -f^c, square = 0. Consequently no folding transition occurs at positive 
K for the square lattice, which remains flat. That Kc is positive for the triangular lattice 
follows directly from the positivity of the entropy at -PC = 0, which makes the triangular 
case more interesting. 
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Fig. 10: Phase diagram of the square lattice folding problem. 

A surprising outcome of our study of the triangular lattice case is the absence of 
intermediate magnetization states between M = and M = 1. Denoting by —Fm{K,0) 
the contribution of the configurations with magnetization M — m to the zero field free 
energy, a sufficient condition for having no intermediate magnetization is that the critical 
field hc,m{K) governing the transition between the M = and an hypothetical M = m 
phase be always larger than hc,i{K) = hc{K), for K < Kc, namely 



with m < 1. The freezing of the system for K > Kc in, say, the M = 1 phase is probably 
a consequence of the intrinsically non-local character of the folding constraint, preventing 
the creation of bounded domains of down spins inside a groundstate of up spins. 

The stability of the M — phase at positive K in the presence of a magnetic field 
h < hc{K) is more surprising. A way of refining the study of the M = phase is to 
introduce a new order parameter, the staggered magnetization 



where the sum alternates between triangles pointing up and down in the lattice. One can 
then distinguish between the disordered folded state Mgt = and a compactly ordered 
folded state Mgt > where the triangles start to pile up. In the presence of a staggered 
magnetic field, we expect the system to behave qualitatively like an antiferromagnetic 



hc,m{K) 



F{K, 0)+Fm{K, 0) 



>h,,,{K) = -F{K,0) + ^K , 



m 
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Fig. 11: Heat capacity ci^\K, 0) for L = 4, 6, 8 and K G [-0.8, 0]. 

(face) triangular Ising model, with a continuous "piling" transition at some negative value 
of the coupling Kc^st < 0. 

Indeed, as explained in 0, the situation is very different when we start to unfold the 
antiferromagnetic M^t = 1 groundstate than when we try to fold the ferromagnetic M = 1 
one. Local deformations of the completely folded groundstate are allowed and enable a 
"low temperature" (large negative K) expansion in terms of a gas of loops of unfolded 
bonds, quite similar to the standard loop gas expansion of the Ising model. As far as 
numerical results are concerned, we only see on Fig. 11 the slow emergence of a peak in the 
specific heat at a value -f^c, st — —0.3. At = 0, an exact solution predicts a disordered 
folded state Mst = with finite staggered susceptibility. 

7. Conclusion 

In this paper, we have derived the phase diagram of the constrained spin system 
describing the folding of the triangular lattice. In the presence of a magnetic field, we found 
a critical line along which a first order transition takes place between a zero magnetization 
phase and the M = 1 pure state, terminating at a triple point {K = Kc,h = 0). This 
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transition persists at zero magnetic field, now driven by the coupling and can be 
interpreted as a first order folding transition between a folded phase and the completely 
flat state of the lattice. The latter is reminiscent of the crumpling transition of tethered 
membranes which is however continuous rather than first order. 

The phase diagrams of Fig. 6 and 10 are very far from the usual unconstrained Ising 
ones. It would be interesting to investigate the role of the local folding constraint by 
applying it gradually to the Ising model, and by looking at the deformation of the phase 
diagram. In this framework, one should be able to follow the evolution of the system from 
a continuous second order phase transition to a first order one. Note that in the square 
case both the Ising and constrained models are special cases of the eight vertex model in 
an electric field, yet to be solved. 
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